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Abstract 

A simulation is performed aiming at checking the 
existence of a well defined stationary state for a 
two dimensional system of driven hard disks when 
energy dissipation takes place at the system bound- 
aries and no bulk impurities are present. 
PACS: 02.70.Ns, 05.60.-k, 47.27.ek 

Bulk dissipation is often assumed to explain sta- 
tionary states in driven systems as in the well 
known example of Drude's theory of electrical con- 
duction where three mechanisms act over a given 
interacting particle system: 

(1) a constant force that accelerates the particles 
in a given direction, 

(2) a thermal bath that should drive the system to 
an equilibrium state absorbing energy excess due 
to the action of the driving and 

(3) an array of bulk impurities that introduce a 
strong chaotic behavior on the particle dynamics. 

The stationary state is characterized by a net 
current of particles in the field direction generated 
by the external work per unit time done by the field 
over the particles equals the heat flux absorbed by 
the thermal bath. The existence of such stationary 
state is physically intuitive: bulk forcing is com- 
pensated by bulk dissipation and it can be seen in 
several computer simulations (see for instance 1). 
Moreover, it is expected that the thermostat model 
used does not influence the system statistical prop- 



erties 011], but some other dynamical properties 
may depend on the particular dissipation scheme 
used [7]. 

Existence of a stationary state is, however, not 
so clear if the action of the thermostat is at the 
system boundaries and no impurities are present: 
the field tends to align the particles trajectories and 
the boundaries introduce a disorder "transversal" 
to the field and this is a bulk versus a surface effect. 

In this similar system studied in hydro- 

dynamics seems to lead to a well defined stationary 
state: a model for the Poiseuille flow in the weak- 
flow regime. There, a group of interacting particles 
are subject to a small external gravitational field 
that drives the particle fiow between two parallel 
plates that are kept at constant temperature while 
strongly interacting via long range forces. Several 
computer simulations of these system confirm that 
the heat generated by a bulk force is efficiently re- 
moved by the thermostats even though the dynam- 
ics in the bulk of the system is conservative: and 
the system evolves tending to a stationary state 
[31 E] . This is different from other studies showing 
thermostats efficiency in cases in which thermostats 
act on the system through its boundaries but the 
driving force also acts on the boundary only, 0]. 

It appears uncertain, more generally, whether the 
boundaries and the likely intrinsic chaotic behav- 
ior of the system would compensate the action of 
a strong external field leading to establish a sta- 
tionary state (the question has been recently again 
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Figure 1: Typical configuration of the model simu- 
lated. The disks in the white part (bulk particles) 
are accelerated in the y direction by a driving field 
E . The disks at the grey boxes act as thermal baths, 
keep constant their overall respective kinetic ener- 
gies (temperatures Ti and T2) for all times. The 
disks interact each other with normal elastic colli- 
sions. The center of the disks may also elastically 
collide with the walls (black lines). 

raised in the literature, [Sj El) and called the 
problem of efficiency of a thermostat mechanism) 
we consider worth studying in other cases the prob- 
lem of whether a thermostat acting only on the 
boundary of the system is efficient enough to ther- 
malize a system subject to bulk driving forces look- 
ing for other similar instances to add to the basic 
result in E] . The latter has been, to our knowl- 
edge, the first to show that such thermostats can 
actually be efficient to remove the heat generated 
by a bulk force even though the dynamics in the 
bulk of the system is conservative: and the system 
evolves tending to a stationary state. 

A thermostat is "efRcient" if it absorbs enough 
energy ( "thermalizes" ) to forbid indefinite energy 
growth of a forced system. 

In this note therefore we consider a system of 
hard disks under the action of a driving field 
and check that, within the range of external force 
strength that we are able to simulate, it appears 
to reach a well defined stationary state. Even for 



"strong" driving fields, although the thermostat 
acts only near the system boundaries, i.e. through 
very short range forces (in fact we consider hard 
core forces) between pairs of particles of the sys- 
tem and of the thermostats and between "mixed 
pairs" of system and thermostats particles. 

The model consists of hard disks confined to a 
unit box and initially placed on a triangular lat- 
tice structure. The disks radius, r, is fixed so 
that the maximum number Nmax, of particles in 

close-packing for a unit surface have a preassigned 

1/2 

value. That is, r = {pcp/NmaxT^) where pcp = 
7r/2V3 ~ .9069 is the close-packing mass density 
for hard disks: we take here N,nax — 10^ and hence 
r = 5.3710-3. 



30 




t 



Figure 2: Evolution of K ^ ^ E^IT" + ^d, 
kinetic energy per particle (top) and particle cur- 
rent J = X^iSi"" ^yi on the y direction (bottom) 
for the bulk disks and for different values of the 
driving field. 

The box will be divided in three parts (see figure 
^ : a central part of width 1 — 2q; ( "bulk" ) with top 
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and bottom identifies ("vertical periodic boundary 
conditions" ) and two equal lateral parts of width a 
("baths"). 

The actual number of disks placed in the bulk 
and in the bath parts are controlled by the cor- 
responding densities of disks: Pbuik and pbath- In 
our simulations we have chosen pbuik = 0.4 and 
Pbath = 0.5. This implies that in our simulations 
the number of disks present in the bulk part is 
Nbuik — 2301, at each bath Nbath — 1318 and the 
total number is then TV = 4937. 

Disks dynamics depends on the sector they are. 
If a disk is in the bulk part it is subject between 
collisions to a uniform acceleration of magnitude 
E in the y direction while in the x direction its 
velocity keeps constant. The disks pertaining to 
the baths move at constant speed along its velocity 
vector. 

In all cases, when the boundaries of two disks 
meet ( "collision" ) they undergo an elastic collision. 
When the center of a disk hits any of the walls 
that define the region in which it is contained, it 
is elastically reflected. In this way we manage to 
keep disks confined at their respective regions and 
particles from the bulk may interact with particles 
in the thermal baths only across the walls. 

Four equidistant walls along the x direction, see 
figure 10, in each bath prevent the disks of the 
bath having a net movement along the y direction 
induced by the interaction with the disks on the 
bulk. 

The disks in each bath keep their total kinetic en- 
ergy constant: K12 — NbathTi^2- This is achieved 
by the following prescription: when a disk from 
the bath collides with another of the bulk, the in- 
crement, A, of kinetic energy (positive or negative) 
that the bath particle suffers is immediately shared 
with the other particles of the bath by rescaling 
their speeds by the factor IS./ Ki^2)~^^'^ respec- 
tively. 

Initially we let the system evolve during 100 7V 
collisions with E = 0: this is empirically sufficient 
to homogenize it spatially; next we turn on the 
driving field. Then, we take measures at intervals 
of N collisions during 10^ iV collisions. We have 
simulated the cases with driving fields E = 0.0001, 
0.0005, 0.001, 0.005, 0.01, 0.05, 0.1, 0.5, 1, 2, 3, 4 
and 5 with Ti = T2 = 1. 

Figure [3 shows typical evolutions of the ki- 
netic energy per particle and the average current 
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Figure 3: Average kinetic energy per particle (K) 
and averaged particle current (J) at the system sta- 
tionary state for the bulk disks. Lines are linear or 
square regression fits of the data. Kj is the internal 
kinetic energy per particle. Inset shows the forward 
discrete derivative of J{E). Error bars are present. 



along the field direction of the bulk disks. After 
some short initial transient, apparently the system 
reaches a stationary state with a well defined cur- 
rent and kinetic energy. 

Figure |31 shows the measured stationary values 
of the kinetic energy per particle and the particle 
current. We see that the hard disk system follow a 
nonlinear current-field response with: J E and 
K oc E"^ for large E while for small fields {E < 1) 
we see Ohm's law with a larger conductivity: J ~ 
1.7 E (see inset in Figure |3J. 

This behavior is consistent with the picture that 
the driving (that tends to align particles) seems 
to generate an intense enough interaction with the 
boundaries that disorders efficiently the bulk par- 
ticles. If we consider the internal kinetic energy, 
Kj — K — J^/2, (total kinetic energy minus the 
kinetic energy of the center of mass in the vertical 
direction) then we find that Ki increases quadrat- 
ically with E exhibiting such disordering effect of 
the boundary interactions which has the effect of 
lowering the conductivity. 

We have also studied the fluctuations distribu- 
tion of J{t) and K{t) around its stationary average 
value. In particular the top of figure ^ shows the 
distribution of the observed values oip{t) — J{t) j J, 
n(p). We see that for £' < 1 the measured distri- 
bution is compatible with a gaussian distribution 
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Figure 4: Top: The measured probability distri- 
bution of p{t) = J{t)/J, n(p), at the stationary 
state for different electric fields. m2{E) is the ob- 
served variance of the p{t) variable. Solid line is 
the normal distribution with zero average and vari- 
ance one. Bottom: Inn(p) vs. (p — l)/m2{E) plot 
for E ^ A, E ^ 5. 



with an average value 1 and a variance 7712 that 
depends on the electric field. However, systematic 
deviations from gaussianity is observed for large 
electric field (see bottom of figure 01). Moreover, 
in figure[Slwe show the variance, m2{E), of J(t)/ J 
and K{t) /K for different values of E. We see while 
7712 for the kinetic energy depends weakly on 
the electric field, the variance for the current decays 
with i? as a power: m2{E) ~ E^^-'^ from E — 0.001 
to E = 1. That is, for i? << 1 a large set of parti- 
cles are able to move in the —E direction. This be- 
havior is strongly suppressed as the field increases 
and, for E > 1 most of the particles move along 
the field. Note that the large error bars due to the 
limited amount of data obtained in this simulation 
obscures the analysis of the large deviations prop- 
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Figure 5: The observed variance m2{E) of J{t)/J 
and K(t)/K versus the electric field. Dashed line 
is guide to the eye showing the power law behavior 
E' 
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erties of p{t) and so we are not able to check the 
fluctuation theorem as proposed in 0. 

Finally we have computed the correlation be- 
tween the current and the energy at the stationary 
state: Cj^k = {J{t)K{t)) / J K - 1. For aU cases 
Cj-K is compatible with zero within error bars. 
That is, on the stationary state both magnitudes 
are uncorrelated and, for instance, the value of the 
system kinetic energy is independent on the sign of 
the current for a given configuration. 

Further investigations could be done by changing 
the temperatures of the thermostats to two differ- 
ent values. 
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